Quaternions in Three Dimensions 
Bob Palais 

Abstract. We construct the quaternion algebra [10] "geometrically" by a three 
dimensional analogue of the classic two dimensional geometric description of the 
complex field. The algebraic description of the multiplication operation in three 
dimensions involves the addition of one term. The construction leads to novel 
methods for implementing and interpolating rotations and understanding their 
topology. 

The bilinear map from ordered pairs of vectors in R 3 to quaternions, T : R 3 x R 3 — > H 
defined by: 

T(v, w) = [v • w, v x w] (1) 

is known to be surjective. (For example, see [3], Theorem 1.) Note that in [3] quaternions 
are viewed as 2 x 2, self-adjoint, complex matrices, while 'vectors' are identified with 
elements of the subspace of such matrices having trace zero. However here will view 
vectors v, w as elements of R 3 and, for the moment, we will pretend that the quaternion 
structure had not been invented and simply view T as a map to Q = [q, q] G R x IR 3 = IR 4 . 
Although T is not injective, if we define an equivalence relation, 

(v, w) ~ (v', w') <^ v • w = v' ■ w' and v x w = v' x w', (2) 

then we get a bijective map T(v, w) = T(v, w) defined on the set G = (R 3 x R 3 )/ ~ of 
its equivalence classes, the inverse images T _1 [g, q] of T. It then becomes natural to ask: 

"Can we introduce a 'multiplication' operation for ordered pairs of vectors that 
makes the non-zero classes, G*, a group, and in such a way that T is a multi- 
plicative homomorphism with some corresponding operation on R — 0? 

Our goal in this note is to show that this question does have an affirmative answer, 
and moreover we will define the multiplication operation by generalizing a very familiar 
geometric construction. When we are done, we will see that its expression as an operation 
on R x R 3 is [g',q'][g,q] ->■ [<?",q"], where 

q" = qq' — q • q' and q" = qq' + q'q Iq'xq. (3) 

This formula is a modern representation of the quaternion multiplication formula that first 
appears in [13]. It is not hard to check that with ei,e2,e3 denoting the standard basis 
for R 3 , if we set 1 = [1,0], i = [0 : ei],j = [0,e 2 ], k = [0,e 3 ], then (3) implies the defining 
relations for quaternions, famously written on a bridge by Hamilton in 1844 [10], 

i 2 = j 2 = k 2 = — 1 and ijk = — 1 

[3]. It follows that ij = k = —ji, jk = i = —kj, ki = j = —ik. Thus our construction 
makes T an isomorphism of G* and H* and gives a purely three-dimensional way to 
interpret (and even 'invent') quaternions. After developing these ideas, we examine in this 
light the implications of viewing known constructions that relate quaternions and rotations. 
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The familiar geometric construction to which we refer combines equivalence classes of or- 
dered pairs of points in Euclidean space to define vector addition, and combines equivalence 
classes of ordered pairs of vectors in the plane to define complex multiplication. (Both of 
these are treated in detail in [1] where they are referred to as 'Chasles' Relation' — though 
[16] questions this terminology. Similar constructions are also used to construct the ring 
of integers from the positive integers and the rationals from the integers.) What we will 
show is that a straightforward generalization of those well-known geometric constructions 
does indeed make T a multiplicative isomorphism. 

Using ||v|| 2 1 1 w 1 1 2 = |v ■ w| 2 + ||v x w|| 2 , if we define ||Q|| 2 = q 2 + ||q|| 2 , we see that the 
equivalence class T _1 [g,q] consists of pairs of vectors orthogonal to q that are separated 
by a fixed oriented angle arccos(j|^yj) and the product of whose lengths is ||Q||. (If q = 0, 
this says v and w are parallel in an arbitrary direction.) An important fact is that any 
two of these equivalence classes 'overlap', by which we mean that there exists a unit vector 
u for which there are representatives of both classes having u as their first element and 
also a unit vector v for which there are representatives having v as second elements. This 
makes it possible, in analogy with vector addition and complex multiplication, to define the 
geometric composition of an ordered pair of our equivalence classes, ((v, w), (v ; , w')) by 
finding representatives such that the second element of the representative of the first class is 
equal to the first element of the representative of the second, and then 'merging' them into 
(v", w") = (v', w') o (v, w), the class of the first element of the first and the second of the 
second. Of course we must show that this geometric merging operation o, is well-defined 
and that it makes the non-zero classes a group isomorphic to the multiplicative group of 
non-zero quaternions. Both of these facts are a consequence of the following theorem. 

Theorem 1. Let v, w, v', w' be vectors in R s with w = v' and ||w|| = 1. Then 
v • w' = (v • w)(v' • w') - (v x w) • (v' x w'), (4a) 

and 

v x w' = (v • w)(v' x w') + (v' ■ w')(v x w) + (v' x w') x (v x w). (46) 

Remarks: Before giving the proof, we note that if we put T(v, w) = [q, q], T(v', w') = 
[g',q'], and T(v, w') = [q",q['\ the conclusion becomes equation (3). Also note that the 
restriction ||w|| = 1 may be removed by an appropriate scaling that puts (4a, 6) in a form 
that highlights their interpretation as a novel vector identity involving just three vectors: 

Corollary. For any three vectors A, B, C G M 3 , 

(A • C)(B • B) = (A • B)(B • C) - (A x B) ■ (B x C), (4a') 

(A x C)(B • B) = (A • B)(B x C) + (B • C)(A x B) + ((B x C) x (A x B)). (Ah') 

While there appears to be an asymmetry, with B plays a distinguished role, A, B, and C 
can be cyclically permuted. We will see that these remarkable vector identities contain 
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the composition laws for (i)the quaternions, H, (ii)the unit quaternions, 577(2), (iii)the 
complex numbers, C, (iv)the plane rotation group 50(2), and (v)the three-dimensional 
rotation group, SO (3)1 

Proof: The identities can be validated by direct brute force comparison of the two sides, 
but since our purpose is to derive them, we prefer a geometric approach. By the orientation 
preserving Givens' QR factorization algorithm for the matrix with columns {v, w, w'}, we 
can always construct a positively oriented orthonormal ordered basis {ei,e 2 ,e 3 } for R 3 
such that ei = w and v is in the e x — e 2 -plane. To motivate (4a, 6)we first assume w' 
is co-planar with v and w, i.e., also in the e x — e 2 -plane. In coordinates with respect to 
this basis, q =< 0, 0, ||q|| >, q' =< 0, 0, |]q'|| > so v =< g, — ||q||, >, w =< 1, 0, > and 
w' =< g', 1 1 q' 1 1,0 >. We confirm by direct calculation that 

q" = v • w' = qq - q • q', q" = v x w' = gq' + g'q (5) 

which we recognize from the geometric meaning of dot and cross products as complex 
multiplication: [cos(0) cos(0') - sin(0) sin(0'), < 0, 0, cos(0) sin(0') + cos(0') sin(0) >]. We 
now let w' be arbitrary. Then, since q' remains orthogonal to w = e 1? it takes the form 
q' =< 0, ||q'|| cos(0), ||q'|| sin(0) > and so w' =< g', | |q'| | sin(^), — | |q'| | cos((p) >. Again 
by direct computation, 

q" = v • w' = gg' - q- q', q" = v x w' = gq' + g'q+ < ||q|| ||q'|| cos(0), 0, > . (6) 

To complete the derivation and proof of (4a, b) we observe that the given data must de- 
termine this cross-product in terms of equivariant quantities and recognize the last term 
of (5) as the non-commutative quantity q' x q. Since in (4a, b) the dot-product is invari- 
ant (Rx • Ry = x • y ) and the cross-product equivariant (Rx x Ry = R(x x y)), they 
remain valid not only when the basis or triple of vectors are simultaneously rotated, but 
also under any rotation of the pairs (v, w) and (v', w') within their respective equivalence 
classes. This completes the derivation and proof of (4a, 6), and guarantees that our geo- 
metric merging operation o is well defined on G, and that with respect to this operation, 
T is a multiplicative isomorphism. QED. 

What we are doing here is viewing quaternion multiplication as the expression of the 
relationships (4a, b) that hold for the dot-product and cross-products of each pair of a 
triple of vectors in R 3 . (Surprisingly, (4a, b) do not appear in standard collections of similar 
vector identities.) While this seems to be novel, if we restrict the vectors in G to lie in any 
two-dimensional subspace, the final term in (46) , (46') vanishes and T becomes equivalent 
to Wessel's construction of an isomorphism between the nonzero elements of (R 2 x R 2 ) / ~ 
and the multiplicative group of complex numbers, C — {0}. We obtain a novel description 
of the three-sphere by restricting the vectors in G to be unit vectors, in which case the 
(B-B) factors in (4a'), (46 ; ) vanish, and T becomes an isomorphism between (S 2 x S 2 ) / ~ 
and the unit quaternions, 577(2) = 5 3 . In this case we may also associate ordered pairs 
of unit vectors with oriented arcs of great circles in the unit sphere 5 2 C R 3 , in which 
case multiplication of unit quaternions may be viewed as 'vector addition' of these arcs. 
The unit quaternion inverse to (v, w) corresponds to (w, v) and the oppositely oriented 
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arc. Algebraically, the dot-product stays the same while the cross-product is multiplied by 
— 1. Under the combination of the above restrictions, the final term in (46), (46') and the 
(B-B) factors in (4a'), (46') vanish, in which case they revert to the standard trigonometric 
addition formulas, and T becomes an isomorphism between (S 1 x S 1 )/ ~ and the unit 
circle group, U(l) 50(2). 

Continuing the analogy with vectors and complex numbers, we note that the additive 
structure of quaternions is also captured geometrically by G, using the fact noted earlier 
that any two equivalence classes overlap with representatives of the corresponding equiv- 
alence classes that have the same first elements: Q = T(v, w),Q' = T(v, w'). Linear 
combinations of the second element are preserved by T because of the bilinearity of both 
dot-products and cross-products, i.e., cQ + c'Q' = T(v,cw + cV). In particular a convex 
combination (and so interpolation) is preserved. But even more is true ! We can show 
that 'SLERP' (quaternion) interpolation [14] on S 3 : 

Q(f) = c(f)Q + c'(t)Q', c(t) = Sin( ^ ( ~ , c'(t) = n = arccos(Q ■ Q') (7) 

can be performed entirely on S 2 . To see this, note that Q(£) = T(v, c(£)w + c'(t)w') 
with the same c(t), c'(t) as in (7). Then since 

Q • Q' = qq' + q • q' = (v • w)(v • w') + (v x w) • (v x w'), 

the Lagrange-Binet-Cauchy identity 

(A x B) • (C x D) = (A ■ C)(B • D) — (A • D)(B • C) (8) 

reduces this to Q • Q' = w ■ w'. In other words c(t) and c'(t) do not change if we replace 
the definition of O in (7) with O = arccos(w • w'). 

This viewpoint also provides insight into the well-known double-covering homomorphism of 
SU(2), the group of unit quaternions, onto the group SO (3) of rotations of M 3 . This takes 
Q G SU(2) to R(Q) G SO (3) defined by conjugation: [c, v] H> Qfc^Q- 1 = [c,R(Q)v], 
or equivalently by the Euler-Rodrigues-Shoemake formula [9]: 

-1*3 

R(Q) = I + 2gJ q + 2J q , where J u := \ u 3 -u x \ . (9) 

Ui 





When we view a unit quaternion as (v, w) G G, this homomorphism take the simple form 
C : (v, w) I—)- C(v, w) = p w Pv, where p u = 2uu T — I is the orthogonal reflection across the 
one-dimensional subspace spanned by the unit vector u. The fact that the double covering 
map C is well-defined on equivalence classes of (2), is a consequence of the relationship 
between rotations R that fix q = v x w, and reflections p u with u • q = 0: 



'Ru 



RpuR" 1 . (10) 
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Then the fact that C is a homomorphism is a consequence of the idempotency of the 
reflections in the unit vector, u, used to compose two classes, (p u ) 2 = I- Because p_ u = p u , 
C(v, w) = C(— v, w) where (— v, w) ~ (v, — w) both represent the negative of (v, w). As 
oriented arcs, they are represented by the arc complementary to that of (v,w), from 
which 2(6 + 7r) = 26 mod 2n also shows that they correspond to the same rotation. The 
remaining fact, that its kernel consists only of the two classes (u, u) and (u, — u) for any 
u £ S 2 follows from the observation that p u u = u and p v v = v and therefore if v is 
not a multiple of u, C(u, v) cannot fix either u or v. Therefore we may view SO(3) as 
equivalence classes of lines though the origin, with a geometric composition analogous to 
o, and an algebraic composition derived from (4a), (46) by identifying classes with opposite 
signs. 

The map C is related to two algorithms for implementing the unique rotation R that takes 
a unit vector ui to another unit vector up ^ — ur and fixes the orthogonal complement 
of their span. The first, motivated by Cartan's 'transvection' isometry [2], has been used 
for an especially efficient implementation of a virtual trackball [5] in the mathematical 
visualization projects, 3D-XplorMath [11] and the Virtual Math Museum [12]. It defines 
s = u; + Uf and with implicit normalization to avoid square roots, then applies R _1 
to the vectors e of the viewing frame in the form p s (p UF e). This highlights the fact 
that the angle between any vector and its image under the rotation C(v, w) is twice 
the angle between v and w. This observation makes it possible to compose two three- 
dimensional rotations "with our fingers" , by composing the corresponding unit quaternions, 
thus forming a virtual SO(4) (or, depending on interpretation, SU(2)) "slide rule". If you 
represent the two planes and half-angles of the rotations using the thumbs and forefingers 
of your two hands, and then find the equivalent configuration so that the interior thumb 
and forefinger overlap, then the exterior thumb and forefinger represent the plane and 
half-angle of the composed rotation! It also explains the difference between the quarter- 
turn rotation about e3 corresponding to the complex i and the half-turn rotation about 
ei corresponding to the quaternion i. The other algorithm explicitly calculates the matrix 
for R in a novel form: 

R = I - — ^— ss T + 2u F u/ T . (11) 
(s-s) 

This can be implemented using only 18 multiplications and one division, and unlike the 
formula in [7], (11) remains valid in any dimension! 

Finally, we observe that the intersection of S 2 with the family of planes in M 3 containing 
the standard basis vector e! whose normals vary at unit speed along the geodesic from e 3 
to e! is 

e(s,t) = sin(t/4)(sin(t/4)ei+cos(t/4)e3)+cos(t/4)(cos(s)(-cos(t/4)ei+sin(t/4)e 3 )+sin(s)e2). 

This explicit one-parameter family of closed paths on S 2 (with base-point ei) starts from 
one turn about the equator in the ei — e2-plane, namely e(s,0) = cos(s)ei + sin(s)e2, 
and shrinks smoothly to the constant path e(s,2-7r) = ei. Then R(s,£) = C(ei, p e (s,t)) 
and Q(s,t) = T(ei,p e ( S] t)) are explicit homotopies of the paths in SO (4) and in S s , 
respectively, between two turns about a fixed axis (e3) and the constant path at the 
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identity. These deformations are the basis of the famous Dirac Belt Trick [15], Feynman 
Wine Glass Trick, and the orientation entanglement relation demonstration described in 
[6] and visualized in [8] and [12]. All of these express, in different ways, the fact that the 
fundamental group of 50(4) is isomorphic to Z 2 , but the current perspective relates them 
to the less mysterious fact that S 2 is simply connected, and also suggests why the same 
construction fails for an odd number of turns or if the axis is constrained to be fixed during 
the deformation, i.e., for plane rotations. 
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